DISCRETE BISPECTRAL DARBOUX TRANSFORMATIONS 
FROM JACOBI OPERATORS 

F. ALBERTO GRUNBAUM AND MILEN YAKIMOV 

Abstract. We construct families of bispectral difference operators of the form 
a(n)T + b(n) + c(n)T~ 1 where T is the shift operator. They are obtained 
as discrete Darboux transformations from appropriate extensions of Jacobi 
operators. We conjecture that along with operators previously constructed by 
Griinbaum, Haine, Horozov, and Iliev they exhaust all bispectral regular (i.e. 
a(n) 7^ 0, c(n) ^ 0, Vn G Z) operators of the form above. 



1. Introduction 

The bispectral problem, as originally formulated by Duistermaat and Griinbaum 
q, asks for a description of all situations where a pair of differential operators in 
the variables x and z have a common eigenfunction ^>(x, z) 

(1.1) L(x,d x )*(x,z) = \(z)*(x,z), 

(1.2) B(z,d z )y(x,z) = 6(x)V(x,z). 

For simplicity we say that L , or B, or 'J, are bispectral when the situation above 
holds. 

The results in || already revealed a number of interesting connections with 
a variety of topics ranging from the Korteweg-deVries equation to the problem 
of isomonodromic deformations for differential operators with rational coefficients. 
Later even more unexpected connections with different areas of pure mathematics 
were found. These include automorphisms and ideal structure of the Weyl algebra 
in one variable g, [|, representations of the Wi+oo algebra g], Calogero-Moser 



system |2J], Huy gens' principle 
(quantized) simple Lie algebras 



traces of intertwiners for representations of 
(the last two in the multivariable case). 
In j(3| all bispectral differential operators L(x,d x ) of second order were classi- 
fied. Notice that if one insists that B(z,d z ) should also be of order two then one 
is necessarily dealing with the Bessel or Airy cases. In that paper very explicit 
use was made of the Darboux transformation mapping a given second order dif- 
ferential operator into another one. When starting from an appropriate bispectral 
L(x, d x ) this was shown to produce another such, with a different B(z, d z ). Wilson 
(23) approached the problem from the viewpoint of commutative algebras of differ- 
ential operators. He classified all maximal bispectral algebras of rank one (which 
by definition is the greatest common divisor of the orders of all operators in the 
algebra). In [l?]] the idea of applying Darboux transformations to commutative 
algebras of differential operators was developed. This allowed for a unification of 
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the apparently unrelated methods in || g3| and an extension of them to the higher 
rank case. Further interesting results in this direction were obtained in |l6| . 

Griinbaum and Haine considered || a discrete-differential version of the above 
problem when the variable x runs over the integer lattice Z and accordingly one 
replaces the differential operator L(x,d x ) by a difference operator 

9 

L{n,T) = Y,h{n)T\ b p (n),b q (n) £ 

i—p 

acting on a function f(n) : Z C by 

(L/)(n) = £>(n)/(n + i). 

i—p 

Following [^2], ^l], we define the support of L(n,T) to be the ordered pair [p,q\. 
Such a difference operator will be called regular if the first and the last coefficients 
b q (n) and b p (n) are nowhere vanishing functions on Z. 

This problem is in a sense a generalization of the problem of classifying orthogo- 
nal polynomials which are eigenfunctions of differential operators. The point is that 
the standard three term recursion relation gives rise to a very special type of differ- 
ence operator, represented by a semiinfinite tridiagonal matrix. In jjj, Griinbaum 
and Haine showed that all instances of difference operators with support [—1,1] 
and second order differential operators satisfying ( |l.l| )-( |L~2| ) result by replacing the 
variable n in the classical cases of the Hermite, Laguerre, Jacobi, and Bessel poly- 
nomials (discovered by Bochner, see Qj) by a variable n + s with n running over the 
integer lattice and e arbitrary (see Sect. 2.2 for precise definitions). The differential 
operator in z is the celebrated hypergeometric second order differential operator 
of Gauss. It is worth noting that the corresponding eigenfunctions ^(n, z) are no 
longer polynomials. For a recent survey of this area, see |D3]. 

Very recently Haine and Iliev classified all maximal bispectral difference algebras 
of rank one E3 . Their result is a beautiful extension of Wilson's work [^3| where 
the Grassmannians associated to Darboux transformations on differential operators 
are substituted with flag varieties coming from such transformations on difference 
operators. Among these some algebras that contain an operator with support [—1, 1] 
were isolated in [p.5| , where they were conjectured to be all of this type. 

The aim of this paper is to make progress in obtaining a discrete-continuous 
analog of the result of |@| , namely a classification of all discrete bispectral operators 
of the form a(n)T + b(n) + c(ri)T _1 (referred to as the extended Bochner-Krall 
problem in jl3| ) . In Jll| the Darboux process was applied to a biinfinite extension 
of the Laguerre difference operators considered in |IJ . A large class of bispectral 
difference operators of the form above was thus constructed and many properties 
of the resulting objects were analyzed in detail. It is fair to say that the results in 



11 1 provide a general treatment of the Laguerre case. The case of Jacobi difference 
operators has so far not been amenable to a similar treatment and only some special 
cases of Darboux maps were proved to preserve the bispectral property. The goal 
of the present paper is to provide such a general treatment in the Jacobi case and 
to state a conjecture for the classification problem above. 
The rest of the introduction describes our results. 
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We take as a starting point the following natural extensions of the Jacobi poly- 
nomials, constructed in : 

p^(n, z) = ( £ + Q: + 1 )» F( _ (n + £ )^ n + £ + a + p + \ ia + l,(l- z)/2). 
\s -\- \ ) n 

Here and later we use F for Gauss' 2-fi hypergeometric function. For negative 
integer values of a, see (2.22). They are no longer polynomials but are still eigen- 
functions of a biinfinite difference operator L Qjj g ;e (n, T) of the form ao(n)T +bo(n) + 
co(n)T _1 and a differential operator B a ^(z,d z ) : 

L a ^ e {n,T)p a /'(n,z) = zp a £ ^(n,z), 

B a> p(z,d z )p^{n,z) = X e (n)p^(n,z). 

The operator L a ^-^{n,T) is obtained by the formal change of variables n i— > n + e 
from the standard (difference) Jacobi operator and is explicitly defined in ( 2.14 ). 
The operator B a p(z,d z ) and the spectral function X £ (n) are given in eqs. (2.10), 

The sets of difference operators that we consider are obtained by the following 
version of the Darboux map starting from the operators L a ^- e (n,T). Let P(n,T) 
be a regular difference operator whose kernel is preserved by L Qj) 9 ;e (n, T). Then 
there exists a (unique) difference operator L(n,T) such that 

(1.3) L(n, T)P{n, T) = P(n, T)L aM , e (n, T) 

which we refer to as a Darboux transformation from L a ^- £ (n, T). The advantage of 
this version is that this L(n,T) is necessarily of the same form as L a ^- £ (n, T), i.e. 
L(n,T) = a(n)T + b(n) + c(n)T _1 for some functions a(n), b(n), c(n), neZ. 

If q(x) denotes the characteristic polynomial of the endomorphism L(n, T) acting 
on the finite dimensional space KerP(n,T), the 

KerP(n, T)cKer 9 (L ai/3 . e (n, T)). 

In view of this it is natural to parametrize the sets of operators L(n, T) by the Grass- 
mannians of special subspaces of Kerq(L a ^- £ (n, T)) that can occur as KerP(n, T). 
Denote the set of difference operators L(n, T) corresponding to characteristic poly- 
nomial q(x) = (x — l) k (x + 1)' by 

ce,/3;e- 

The operators in D^' E are the main objects of study in this paper. Their explicit 
form is given in Sect. 3.3. Restricting to q(x) with roots at ±1 guarantees that 
L(n, T) will have rational coefficients. This is an important feature of bispectral 
operators. See || in the differential case. 

It is an easy consequence of (1.3) that the function 

(1.4) *(n, z) = P(n, T)p^{n, z) 

is an cigenfunction of the operator L(n, T), namely we have 

L(n, T)*(n, z) = z^(n, z). 

Our main result is: 

Theorem 1.1. The difference operators L(n,T) from the sets "D^'ji-s are bispec- 
tral [or more precisely the functions ^(n, z) ( |l.4| ) are eigenfunctions of differential 
operators in the variable z) in the following cases 
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1. ) a € Z and k < \a\, I = 0, 

2. ) P G Z and Z < \/3\, k = 0, 

3. ) a, f3 eZ andk< \a\, I < \/3\. 

When e = 0, and k = 1 and (or) Z = 1 these results were established in [|l8[ 
All this work starts with the classical paper of H. L. Krall, |19[. 



The proof of Theorem 1.1 is based on a general result of Bakalov, Horozov, 
and Yakimov |j| which guarantees that a Darboux transformation preserves the 
bispectral property under some conditions on the operator P(n,T). We will soon 
see that its application to the present situation is highly nontrivial and requires, in 
particular, an intrinsic characterization of a space of difference operators. 

We will need some notation from ||, see Sect. 4.1 for more details. Denote 
by Ba,(3-,e the algebra of difference operators S(n, T) with rational coefficients for 
which there exists a differential operator G(z, d z ) (also having rational coefficients) 
satisfying 

(1.5) S(n, T)p«<P(n, z) = G(z, d z )p^(n, z). 

All such operators S(z, d z ) form a "dual" algebra B' a g. e . The map 

l>:B,,: b(R(n,T))=S{z,d z ) 

is an antiisomorphism of associative algebras. Let /C a ,/3 ; e and K' a o. e be the subalge- 
bras of B a ^- E and B' a ^ consisting of rational functions. Bispectrality oip"'^(n, z) 
is equivalent to /C a ,/3 ; e and K! a g. s being both nontrivial. Finally we arrive at the 
most important object for our consideration, namely the space 

TZa,/3;E = {(v( n )y 1 Po(n,T) | fjb(n) € lC a ^- s , P (n,T) £ S Qj/3;e ,and the operator 

(/j,(n)) _1 Po(n, T) does not have poles atn 6 Z}. 

According to Theorem 1.2 of ^(n, z) is an eigenfunction of a differential operator 
in the variable z, if 

P{n,T)en a , P;s . 

Thus to prove Theorem 1.1 we need a good description of the space K. a ^ ;s which 



can be used to check whether the operators P(n,T) from (|L3|) belong to TZ a ,/3-e- 
This is the hardest step in our paper. Let A denote the algebra of abstract difference 
operators with rational coefficients of the form YH= P bi(ri)T l with rational functions 
bi(n) (possibly having poles in Z). The key point of our approach is to consider the 
involution / of A acting on rational functions h(n) by 

(Ih)(n) := h{-{n + 2s + a + /3 + 1)) 

and on the shift operator T by I(T) := T _1 . In Sect. 4.2 we prove that H a ,p-e 
consists of those difference operators from A that do not have poles in Z and after 
conjugation with the function 

(e + a + l) n 
^ = (* + !)„ 

become /-invariant. 

The final step of the proof of Theorem is to show that the hypothesis guaran- 
tee that the kernel of the operator P(n, T) (defining L(n, T)) has a basis of functions 
f(n) for which the ratio f(n)/(j)(n) is an (almost) /-invariant rational function in 
n. This is done in Sect. 5.1. Finally Sect. 5.2 recapitulates the strategy of the proof 



DISCRETE BISPECTRAL DARBOUX TRANSFORMATIONS 



5 



of Theorem LI for the special case of the set T>^'^ e . The reader may find it useful 



to consult this section while reading the paper. 

Let us also note that in the case of Laguerre polynomials the situation simplifies 



a lot due to a presense of a relation of the type (1.5) with a difference operator 
S(n,T) of the form s\{n)T + sq(u) and a first order differential operator G(z,d z ) 
(see expressions (2.3) and (2.8) in |ll]]). It is not hard to show that as a consequence 
of this the analog of TZa.p-e m that case is simply the space of difference operators 
with rational coefficients. 

Comparing with the differential case it is natural to conjecture that all sec- 
ond order regular bispectral difference operators {i.e. having support [—1,1]) are 
exhausted by the families of operators constructed in M, IllL [15j and in this article. 



The operators in |15 are obtained as Darboux transformations from the operators 
ia.ftei",! 1 ) for half integer values of the parameters a, (3 and are the analogs of 
"KdV family" in the differential case ||. 

For later use we introduce some convenient notation. If f{n) : 1 — > C is a nowhere 
vanishing function and D\(n,T), D 2 {n,T) are difference operators we denote 

Ad /(n) Z?i(n,T) := f{n)D 1 {n,T)f{n)- 1 , 

ad D2( „, T)J Di(n, T) := D 2 {n, T)D x {n, T) - D x {n, T)D 2 {n, T). 



2. BlINFINITE JACOBI OPERATORS 

In the first part of this section we review some properties of the classical Jacobi 
polynomials p"'^(z). The second one discusses certain functions p"'^(n, z) which are 
eigenfunctions of biinfinite analogs L a ^. e (n,T) of the Jacobi difference operators. 
The third part describes Darboux maps between the operators L at p. £ (n,T) with 
shifted indices a, (3. 



2.1. Jacobi polynomials. The Jacobi polynomials are the orthogonal polynomi- 
als for the measure (l~z) a (l+zfd z on the interval [—1, 1], (a, (3 > —1), normalized 

by 



^(i),r»(" T "UEV 



n + a 
n 

They are given by 

(2.6) V a J{z) = (^ + n F(-n, n + a + /3 + l;a + l;(l- z)/2) 

where F (a, b; c; x) denotes the Gauss' hypergeometric function. The reader can 
consult ||^, pp. 209-217] for other explicit formulas and a list of major relations 
for p^P{z). Let 



(2.7) p a >f\n,z) 



\p%' p {z), for n S Z> 
| 0, for n £ Z<o 



Now p a, @(n, z) are functions of a discrete parameter n and a continuous parameter 
z. They satisfy a three term recursion relation 

(2.8) L a ^(n, T)p a ^(n, z) = zp a ^(n, z) 



(i 
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where L a ^(n,T) are the difference operators 
(2.9) 

T , T , 2 (n+l)(n + q + /3+l) , /3 2 - a 2 



(2n + a + /? + l)(2n + a + /3 + 2) (2n + a + 0){2n + a + + 2) 

2(n + a)(n + /3) r _ t 
(2n + a + 0)(2n + a + + 1) 

called Jacobi operators. In addition p a '^(n, z) are eigenf unctions of the differential 
operators B a> p(z,d z ) given by 

(2.10) B a ,p(z, 8 Z ) = (z 2 - 1)9 2 + {a-f3+{a + {3 + 2)z)d x , 
i.e. 

(2.11) B a , p (z, d z )p a ^(n, z) = \{n)p a ^{n, z) 
for 

(2.12) \{n) = n(n + a + (3 + 1). 



In view of ( pTq ) and (2.11), p a >P(n, z) are discrete-continuous bispectral functions 
and L a ^(n,T), B(z,d z ) bispectral difference (differential) operators. 

2.2. The functions p"'^{n, z). In a study of the relation between the so called 
"associated Jacobi polynomials" and the discrete-continuous bispectral problem 
Griinbaum and Haine, see @, [h], |l2), introduced the functions 

(2.13) p»>P(n, z) = ^±£±^ F(-(n + e),n + e + a + + 1; a + 1; (1 - z)/2) 

(£ T" l)n 

(n G Z, z G C, |z| < 1) defined for those e, a, (3 G C such that a ^ Z<o, and e ^ Z<o, 
e + a Z>o- We will see later that the first restriction can be lifted. 

The functions p"'^(n,z) are no longer polynomials but satisfy relations, similar 
to the ones for p a 'P(n 7 z). In particular they are eigenfunctions of the following 
difference operators with support [—1,1] 

, 914 s r f T x_ 2(n + s + l)(n + e + a + l3 + l) 

( ' h »,0;e{n, (2n + 2e + a + /3 + 1) (2n + 2e + o; + /3 + 2) 

/? 2 - a 2 



(2n + 2e + a + /3)(2n + 2e + a +/3 + 2) 
2(n + e + a)(n + e + /3) 



(2n + 2e + a + /3) (2n + 2e + a + (3 + 1) " 

and of the differential operators B a ^(z, d z ), eq. ( 2.10 ). The corresponding relations 
are 

(2.15) L a ^ e (n, T)p?f>(n, z) = zp^(n, z), 

(2.16) B a ,p(z, d z )p^(n, z) = X e (n)p^(n, z). 
where 

(2.17) X E {n) = (n + e)(n + e + a + f3+l). 
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The difference operators L a ^ ;e (n, T) will still be called Jacobi operators. Further 
we will only deal with the case when they are regular, i.e. when their coefficients 
of T and T _1 do not vanish for n G Z. This amounts to the conditions 



Z. 



(2.18) e,e + a,e + P, e + a + f3, 2e + a + /3 • 

It may be useful to stress here that these will eventually be the only restrictions on 
our parameters a, /?, e. 

The operators L a ^. £ (n,T) do satisfy certain "transformation properties". For 
instance the following relations hold 

(2.19) L-a,-/3,s+ a +p(n, T) = L a> p ie (n, T), 

(2.20) Ad(_ 1 )ni (3)a!i£ (n, T) = Ad(_ 1 )nZ_ ( 3 ) _ a , iE+Q , +/3 (n ) T) = -L a .p. £ (n, T). 

It is tempting to use ( 2.19| ) to limit attention to the case a > 0. However, this 
would eventually bring an undesirable degree of asymmetry in the treatment of 
the parameters a and /3. For this reason we prefer to introduce the appropriate 
functions p" ,/3 (n, z) for a e Z <0 (and e, e + /3, e + a, e + a + (3 ^ Z) by using ( [2.13 ) 
and recalling, see [^0[ p. 38] that for m e Z>o 

(2.21) 

1 „, , ^ (a)m+l(&)m+l 



lim 



r(c) 



F(a, b; c; z) 



(m+ 1)! 



-z m+1 F(a + m+ l,b + m + l;m + 2;z). 



We see below that this leads to the following expression for p"'' (n, z) with a £ Z 
(as long as ( 2.18j ) is satisfied) 

(2.22) 

(e + /?+!)„ (1-^)-° 



<o 



C- 



F(-(n + e + a),n + e + + 1; -a + 1; (1 - z)/2) 



{e + a + (3+l) n 2- 
where the constant C = C(a,(3,s) is explicitly given by 

(-1) Q (-e)_ a (e + a + P + l)- a 



C = C(a,/3,e) = 



(-a -1)1 



(-a)! 



It is easy to check that the assumptions (2.18) imply that C(a, /3,e) is well defined 
and does not vanish. 

The expression above can be derived by a continuity argument using (2.13) and 
(2.21) when a approaches a value in Z<o- To see this it is important to notice that 
for a 6 Z<o the identities 

(s + a+l) n (-£)-« . , (e + /?+!)„ (n + e + a + (3 + 1)_ Q 



and 



(e + a + /3 + 1) T 



(£+l)n (-(n + e))- 
allow one to rewrite the factor 

(e + a + l)„ (~(n + e))- a (n + e + a + P - 



(e + a + f3+ l)_ c 

•1)-Q 



(2.23) 



(e + l)n 
(-a 



D! 



(-a)! 

(e + /3 + l) r 



(-1)" v ' (e + a + /3+l)„ 

which except for the first constant is the factor in front of ( [2.22D , 

Then conditions ( ^.18 ) guarantee that p"'^(n,z) is well defined (see (2.13) and 
(fT2§)) and satisfie s Wl§ an d <\2 16| ). It was proved in |9) that the space of com- 
mon solutions of (2.15) and ( |2.16 ) in a domain f2cC, not containing ±1, is two 
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dimensional. Notice also that ( [2.18 ) excludes, in particular, the original operators 
L a ,i3{n, z) (e = 0) since their leading coefficient vanishes for n = — 1. 

Finally we explain how (2.15) follows from fl2.8p . Conjugate the operator L a ^(n, T) 
with ( n + a ) = The resulting differ ence operator has rational coefficients 

and the eigenfunction F(—n,n + a + j3 + 1; a + 1; (1 — z)/2), cf. (2.6). The op- 
erator obtained from it by the formal change n i— > n + s has the eigenfunction 
F(—n — e, n + e + a + (3 + 1; a + 1; (1 — z)/2) and all we need to do is conjugate 
it with (e + l) n /(e + a + 1)„. The result is the operator L a ^ ;e (n,T) which proves 
(fH|) in the case a ^ Z<o- The case a £ Z<o follows from the definition ( |2.13D 
using the limit ( [2.21 ). 

2.3. Darboux maps between Jacobi operators. There are four difference rela- 
tions connecting the values of the Jacobi polynomials p Q:/3 (n, z) with shifted indices: 

\2n + a + p 2n + a + [3 J 



p^(n,z) = -^ I ( 



2(n+ 1) 



-T — 



2{n + a) 



and 



a,/3-l 



(n, z) 



2n + a + + 1 2n + a + (3 + 1 
n + a + (3 n + < 



(n,z), 



T- 1 p a >P(n,z), 



p a >f\n,z) 



1 



2n + a + f3 2n + a + (3 

2(71 + 1) _ 2(n + f3) 



-T 



z + 1 \2n + a + (3 + 1 2n + a + /3 + l 



-,a,/3-l 



(n,z), 



(see for instance, [£C| eqs. pp. 209-219]). Similarly to the proof of (2.15) at the 
end of the previous subsection, one shows the following analogs of these identities 
for p®'P(n, z) 



(2.24) 

pf-Wfaz) = Dl(n,T)p°-P(n,z), p° +1 ' (n,z) 
(2.25) 



2-1 



D«(n,T)p^(n,z), 



= D^(n,T)p^(n,«), p^ +1 (n,z) = J—D f l(n,T)p^(n,z), 



z + 1 



where the operators D±(n,T) and D±(n,T) are given by 



D°L{n,T) 
D+(n,T) 
D P _{n,T) 



Di(n,T) 



e + a 
a 

a + 1 



n + e + a + (3 



n + e + (3 



2n + 2e + a + (3 2n + 2e + a + (3 

2{n + e + l) 2(n + e + a + l) 



T 



e + a + lj \2n + 2e + a + (3 + 2 2n + 2e + a + / 3 + 2/' 
/ n + s + a + (3 



+ 



\2n + 2e + a + f3 2n + 2e + a + (3 



-T~ 



2(n + e + P + l) 
2n + 2e + a + (3 + 2~ ' 2n + 2e + a + f3 + 2 



2(n- 



1) 



-T + 



The constant in (2.22) was chosen to make the relations ( 2.24 )-( 2.25) hold for all 
a E C. We show only the dependence on the index a of the operators D±{n,T) 
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because the index [3 is unchanged in both sides of eqs. (2.24), similarly for the 
operators D±(n, T). Eqs. ( 2.24 )— ( 2.25| ) and ( |2.15| ) imply the following factorizations 



(2.26) L a ,p, e (n,T) - 1 = D^n, T)D«(n, T) = D" +1 (n, T)D"(n, T), 

(2.27) L a ^ e (n,T) + 1 = D^ 1 (n,T)D _(n,T) = D f ^ +1 {n,T)D f l(n,T). 

Hence the operators L a ±i ^; e (n , T), L a ^±i- £ (n,T) are Darboux transformations 
from L a .p- £ (n, T) and eqs. ( 2.24 ), ( |2.25| ) represent the Darboux maps p £ '@(n, z) \— > 



p^ hfi (n, z) and p%> (n, z) i-> p^ 1 ^, z). 



3. Darboux transformations from Jacobi operators 

The first part of this section contains some general facts about discrete Darboux 
transformations in the form in which they will be used later (see, for instance, [ p5| 
for the differential case) . The goal of the second part is an explicit description of the 
kernels of the operators {L a ,p;e — ^) k {L a ,p;e + 1) < Based on it, in the third part we 
construct Darboux transformations from L a p- e {n,T) which are the main objects 
of study in the rest of the paper. The conditions ( 2 . 1 8| ) are assumed throughout 
Sects. 3.2-3.3. 

3.1. General remarks on Darboux transformations. One says that the differ- 
ence operator L(n, T) is obtained by a Darboux transformation from the difference 
operator Lo(n, T) if there exists an operator P(n,T) such that 

(3.1) L(n, T)P{n, T) = P(n, T)L {n, T). 
Assume that Lq(u,T) has an eigenfunction ^o(n,z), i.e. 

(3.2) Lo(n, T)* (n, T) = flb(z)*o(n) 
for some function go(z). Then 

V(n,z) := P(n,T)t> (n,z) 
is an eigenfunction of L(n, T) : 

(3.3) L(n,T)tf(n,T) = so(z)¥(n). 

The map *f?o(n,T) i— > *&(n,T) is also called a Darboux transformation. 

An important feature of the transformation (3.1) for a regular difference operator 
P(n,T) is that the operator L(n,T) has the same support as Lo(n, T). Besides this 
L(n,T) is regular if and only if Lo(n,T) is regular. 



Given a difference operator L$(n,T), all transformations of the type (3.1) with 
a regular difference operator P(n, T) can be described in terms of the kernel of 
P(n,T). 

Proposition 3.1. (i) For a regular difference operator P(n,T) there exists an op- 
erator L(n,T) for which (3.1) holds if and only if 

(3.4) L (n, T)(KerP(n, T))cKerP(n, T). 



The operator L(n,T) satisfying (3.1) is unique 



(ii) Let P(n,T) be a regular difference operator satisfying (3.4) and q(x) be the 
characteristic polynomial of the linear map Lo(n, T) acting in the space KerP(n, T). 



10 



F. A. GRUNBAUM AND M. YAKIMOV 



Then KerP(n,T)Cq(Lo(n,T)) and there exists an operator Q(n,T) such that 

(3.5) q(L (n,T)) = Q(n,T)P(n,T), 

(3.6) q(L{n,T)) = P(n,T)Q(n,T), 

Note that the kernel of a regular difference operator P(n, T) is finite dimen- 
sional. More precisely, if P(n,T) has support [7711,7712] for some mi G Z, then 
dimKerP(n, T) = rri2 — mi. For any j G Z the map 

(3-7) / -» (f(j + 1), . . . , f(j + 7772 - 777!)), for / I Z -» C 



provides an isomorphism between KerP(n,T) and C m2— mi . 

The transformation Q(?7, T)P{n 1 T) 1— > P(n, T)Q(n, T) is a more traditional ver- 
sion of the Darboux map. Although it is a special case of the transformation 
Lq(u,T) 1— ► L(n,T) from eq. ( |3.l|) , Proposition |3.l| shows that there always exists 
a polynomial q(x) for which q(Lo(n,T)) 1— > q(L(n,T)) is a Darboux map in this 
sense. 

Proo/ 0/ Proposition pp. (i) If P(n,T), L(n, T) satisfy @) and / (77) e KerP(77,T) 
then 

P(n,T)(L (7i,P)/(77)) = L(n,T)P(n,T)f(n) = 
which proves ( |3.4[ ). 

In the opposite direction, let us notice that a comparison of the coefficients of the 
two sides of eq. (3.1) for a fixed value of n gives a finite system for the corresponding 
coefficients of the unknown operator L(n, T) having the same support as Lq(ji, T). 
One shows that it has a solution using the standard linear algebra fact that for a 
finite matrix A the system Au — b has a solution if and only if v l b — 0, V?j 6 KerA f . 
In the particular case which we consider the last condition is fulfilled because of 

(0)- 

The regularity of the difference operator P(tt,, T) implies the uniqueness of the 
operator L(n,T) satisfying (3.1). Indeed if there are two such operators L(n,T) 
and L'{n,T) one can subtract the resulting equalities ([O]). This gives (L(n,T) — 
L'(n, T))P(n, T) = which is a contradiction. 

(ii) The relation KerP(n,T)cq(Lo(n,T)) follows from the definition of q(x). 
Similarly to part (i), this implies the existence of an operator Q(n,T) satisfying 
Eqs. Q and (|J) imply 

q(L(n,T))P(n,T) = P(n,T)q(L (n,T)) = (P(n,T)Q{n,T))P(n,T) 



and as a consequence of this (3.6). □ 
A regular difference operator is reconstructed from its kernel by the following 
lemma. 

'erence operator with support 
Span{/W(n)}2. 1 where m = 



Lemma 3.2. Assume that P(n,T) is a regular 
[7711,777,2] and leading coefficient 1. Let KerP(n, T) 
7772 — TTi2 ■ Then the function 

det(n) := det(f w (n-j)) 

does not vanish for n £ Z and 

1 /Wfn + TTH) ■■• /< 

(3.8) P(n,T) 



i,j—l : mi 



det(n) 



/(i)( r 



777 2 ) 



(n 

/M( n 



mi) 



m 2 ) T m2 
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where the determinant is expanded from left to right ( the shift operator T does not 
commute with function multiplication). 



Proof. The fact that the map ( p.7|) is an isomorphism between KerP(n, T) and C m 
implies that det(n) does not vanish for n G Z. Clearly the functions f^'(n) belong 



to the kernel of the operator in the r.h.s. of (3.£). It has leading term 1 and the 



nonvanishing of det(n) implies (|3.8|). □ 



Remark 3.3. The composition of two Darboux transformations Lo(n, T) t— > L\(n,T) 
and Li(n,T) i— » Li(n^T) of the type (3.1) is Darboux transformation Lq(ji,T) \— * 
L2(n,T) of the same type. Indeed if 

L % {n,T)P % (n,T) = P i (n,T)L i _ l (n,T), i= 1,2, 

then 

L 2 {n,T)P 2 {n,T)P 1 {n,T) = P 2 {n,T)P 1 {n,T)L Q {n,T). 

3.2. Description of Ker(L Qj/ 3 ;e — l) k (L a ^. e + l) 1 . The main idea is to first find 
some functions (p(n, z) (depending on a, (3, and e) such that 

(3.9) L a ,/3;e(n, T)(p(n, z) = zip(n, z) 
and then to consider the derivatives 

^£(n)^^dMn,z) 
They satisfy 

(3.10) {L a ^ e {n,T)Tl)^{n) = <pi~ l \n), VieZ>„ 
with ip± ^ (n) — 0. As a consequence of this 

(L«,/j !e (n, T) t 1) V± -) (n) = 0, V i G Z >0 , j = 0, . . . , * - 1. 

Before stating the results from this subsection we recall a relation for the hy- 
pergeometric function that is a consequence of Gauss' relations between contiguous 
hypergeometric functions. Denote F — F(a,b; c; (1 — z)/2), TF = F(a — 1,6 + 
l;c; (1 - z)/2), and T~ 1 F = F(a+l,b- l;c; (1 - z)/2). Then for c ^ Z< 

mil 2(c-a)b 2( a + b-l)(-2 C + a + b+l) 

1 ' J (&-a)(&-a + l) (6-o-l)(6-o+l) 

-T~ l F = zF. 



, i G Z> . 

z=±l 



(6-a)(6-a - 1)" 

This can also be checked directly using the standard expansion of F(a, b; c, x) for 
\x\ < 1, c ^ Z< 

(3.12) F{a ^ c . x) = f^(4A x j. 
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Lemma 3.4. The four functions 
(3.13) 

(e + a + 1) 



(3.14) 

i/j+(n,z) 
(3.15) 

ip-(n,z) 
(3.16) 

ip-(n,z) 

satisfy 
(3.17) 



F(-(n + e),n + E + a + /3 + l;a + l; (1 - z)/2), 

F(-(n + e + a + /?), n + e + 1; -a + 1; (1 - z)/2), 
F(-(n + e), n + e + a + j3 + l;/3 + 1; (1 + z)/2), 
F(-(n + e + a + /?), n + e + 1; -(3 + 1; (1 + z)/2) 



(e + 1)„ 

(e + /3+ l) n 
(e + a + /3 + 1) 

(-!)"(£ + /? + 1) 
(e + 1)„ 

(-!)"(£ + a + l) 
(e + « + /?+ 1)„ 

(L a ,0:e(n,T) ~ z)tp ± (n,z) = (L a ^. E (n,T) - z)ip±(n,z) = 0, 



provided that a £ Z<o (Z>o) for ip+(n, z) (ip + (n, z)) and (3 ^ Z<o (Z>o) for tp-(n, z) 
(i)-(n,z)). 



Note that the assumptions (2.18) guarantee that the denominators of the first 
factors of tp±(n, z) and if>±(n, z) do not vanish. 

Proof. The relation ( p. 17 ) for (p + (n,z) holds because <p+(n,z) = p"'@(n, z). To 
check the one for ip+(n, z), we conju gate L at p ]e (n, T) by (s+/3+l) n /(e + a+f3+l) n 
(the factor in front of the r.h.s. of ( 3.14 )). 
The result is 



^(e+P+l) n /(e+a+l3+l) n L a ^ ]£ (n,T) 

_ 2(n + e + \){n + e + (3 + 1) 

~ (2n + 2e + a + [3 + l)(2n + 2e + a + [3 + 2) 

(3 2 -a 2 



T 



(2n + 2e + a + (3){2n + 2e + a + (3 + 2) 

2(n + e + a)(n + e + a + (3) 
(2n + 2e + a + (3){2n + 2e + a + /3 + 1) 



This is the difference operator from the l.h.s. of (3.11) with a = — (n + e + a + /?), 
b = n + e + 1, and c = —a + 1 which gives the proof of (3.17) for ip+(n, z). The 
cases of tp-(n, z) and ip-( n -i z ) are handled in a similar fashion. □ 

Next we consider the derivatives of ip+(n, z), ip+{n, z) at z = 1 and of f-(n, z), 
%l>- (n, z) at z = —1 : 



^(n) := ^a*V±(",^) 
z! 



z=±l 



z=±l 



i £ Z>o (with the restrictions on a and (3 made at the end of Lemma 3.4). 
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Using the expansion ( 3.12| ) of the hypergeometric function, we obtain the follow- 



ing explicit formulas for (p± in) and (n 



(3.18) 


(p\'{n 


(3.19) 




(3.20) 


(i)/ 


(3.21) 




We define 





(e + a + l) n (-(n + e))j(n + e + a + fi + l)i 

(e + l)n 
(e + /3 + 1)„ 



(-2)^!(a+l)i 
(-(n + e + a + P))i{n + e+ 1), 





+ a + 


/? + 




(-2)<i!(-a + l)< 




( 


1)»( £ 


+ /3 


fl)« (-(n 


+ e))i(n + e + a + /3H 


-1). 










2*i!(/8 + l) i 




(- 


l)"(e 


+ a 


(-(n 


+ e + a + /3)) l (n + e- 





(e + « + /? + l) n 



2*t!(-/3+l)j 



V>+(n)) for a G Z <0 (a £ Z> ), i < |a| by jSTl^) , ( ^19| ) and 
<p®(n) (V'-(n)) for /? G Z <0 (/? G Z> ), i < |)8| by ( p(i| ), ( gjg ). (Note that these 
cases were excluded in Lemma 3.4.) 

Theorem 3.5. Assuming ( |2.18[ ) i/ie following relations 



{L a> p. e {n,T)-l)^\n) = ^>{n), 
(i a , /3;e (n,T)-l)^i i) (n)=^" 1) (n), 



(3.22) 
(3.23) 

/ioZ<i /or a// i G Z>o if a and for i = 0, . . . , |a| — 1 if a E Z. Similarly one has 

(3.24) (i a , te (n,T) + l)^(n) = ^V), 

(3.25) (L a!/3;e (n,T) + l)^ W (n) = V-^V), 

/or aZZ i G Z> i/ /? ^ Z and /or i = 0, . . . , - 1 i/ /3 6 Z. (We sef (n) = 
V4 _1) (") = 0.) 

XTie kernels of (L a ^. E (n,T) — l) k and (L a ^- £ (n,T) + 1)' are given by 

(3.26) Ker(L Q , ftE (n,T) - l) fc = Span{^(n), ^(n)}^ 1 , 

(3.27) Ker(i Qi/3;e (n,T) + 1)' = Span^V^-V)}^ 

for k<\a\ if a G Z, /or Z < if (3 £ Z, and /or aZZ fe, Z > if a, /3 £ Z. 

Proof. In the case a ^ Z, the functions ip + (n, z) and ?p + (n,z) are well defined. 
From the remark in the beginning of this subsection it follows that ( 3.17 ) and the 

definitions of <p+\n), ip^ (n) imply ( 3.22 ) , ( 3.23[) . The case a £ Z, i < \a\ follows 

by continuity on a. 

The inclusion D in ( 3.26 ), ( 3.27 ) clearly follows from ( 3.22] )-( 3.25 ). Because 
(L ai p-e{n,T) — 1) is a regular difference operator with support [— k,k], to prove 

( |3.26| ) it suffices to show that the functions y>+ (ra), ip+(ri), i = 0, . . . , fc — 1 arc 
linearly independent. 
Assume that 

feo 

(<W+ (n) + (n)J = 0, Vn G Z 

i=0 

for some complex numbers do, . . . , afc , 6o, . . . , Ofc , such that eife ^ or &fc ^ 
(fc < fc - 1). Applying {L a ^. e {n, T))^' 1 to this equality and using ( |3.22| ), ( |3.23[ ), 
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we get 



ak p+ (n) + bk ip+ } (n) = 0, Vn e Z, 



(3.28) 



( e + q + l) n 
(e + l)n 



-&fe 



(e + /3+l)n 
(e + a + /3 + 1), 



Vn G Z. 



For n = this gives a,k = — &fc„ 0). Dividing the two sides of eq. fl3.28j ) for two 
consecutive values of n, we get 

(e + a + n)(e + a + (3 + n) = (e + n)(e + + ri), V?i G Z 

This gives a = which is a contradiction. Eq. ( 3.27 ) is proved analogously. □ 

Remark 3.6. It is clear that 

Ker(L ai/9 . 8 - l) fe n Kev(L a ^ £ + I) 1 = 0. 

Therefore 

Ker(L Qi/3;e - l) fc (L Q , ft£ + 1)' = Ker(L Qi/3;e - l) fc © Ker(L Q , fte + 1)' 

and Theorem [D] describes the kernel of the operator (L a ^- e — l) k (L a : p ;e + l) 1 in 
the cases specified there. 

3.3. The sets T>^i: e of Darboux transformations from L a ^. E (n, z). Let us 
fix two nonnegative integers k and I and choose 2(k + I) complex numbers 

A.\ , i ^ ~~ 0, . . . . k lj 
C h D h j = 0,... ,1-1. 
If k > (I > 0) we will assume a ^ -fc + 1, . . . , k - 1 (/3 ^ -Z + 1, . . . , I - 1). Set 
'Ei.=o(A-vJ" r) (n)+^i <_r) (n)), fori = 0,... ,fc-l 



E;=o(C^r K - r '(n)+D r ^- K - r \n)), for i = k, . . . , k + I - 1 
The values of the parameters A,B,C,D G C for which 

(3.29) det(n) = det(/»(n + i))-,+Lo,-fc-i ^ 0, Vn G Z, 
will be called admissible. For such values we define the operator 

1 /(°)(n-fe-0 ••• /(*+'-D(n-fc-J) T-( fc +') 

(3.30) P(n,T) 



det(n) 



/(°)(n) ••■ /(^-^(n) 1 

By expanding ( 3.30|) along the last column one sees that the term of T~( k+l ^ is 
given by 

det(n+l) 
det(n) 

hence P(n, T) is a regular difference operator. As a consequence of properties 
(p^)-(|3!25|) we obtain 

(3.31) (L a , l3;e (n,T)-l)fM(n)=0, (L a , ff . E (n, T) - = /^(n) 
for i = 1 , . . . , fe — 1 and 

(3.32) (L^(n, T) + 1)/W(n) = 0, (i ai/J . B (n, T) + !)/« - (n) 
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for j = k + 1,... ,Jfe + I - 1. Thus KerP(n,T) = Spanf/^n)}^ -1 is preserved 
by L a ^. £ (n,T) and according to Proposition |3.l| there exists a difference operator 
L(n,T) with support [—1, 1] such that 



(3.33) 



L(n,T)P(n,T) = P{n,T)L a , p . c (n,T). 



The set of all difference operators L(n,T) for admissible values of the parameters 
A, B, C, D will be denoted by 



V 



(k,l) 



All operators L(n,T) £ ^a 'a-e are DaAoux transformations from L a ^. £ (n,T) 
and k,l refer to the multiplicity of the eigenvalues 1 and —1 of L a p- e (rt,T) in 
KerP(n,T), see eqs. (3.31) and (3.32). (Recall from part (i) of Proposition 3.1 
that L a ,/3; £ (n,T) preserves KerP(n,T)). Every L(n,T) £ D^*o. e is a regular differ- 
ence operator with eigenfunction 



(3.34) 

more precisely: 
(3.35) 



*(n,z) = P(n,T)p«<P(n,z), 



L(n, T)^(n, z) = z^(n, z). 



The admissibility condition ( 3.29| ) holds for almost all values of A, B,C,D 6 Z. 
The complement of the corresponding set in C 2( - fc+ ^ consists of the zer os of count- 
ably many polynomials, obtained from det(ro) for fixed n € Z (recall ( |3.29 )). The 
latter do not vanish identically due to the linear independence of the set of func- 
tions {tp± (n)}^'^ U {ip±^ (n)} l ~} (see the proof of Theorem 3^5) and the regularity 
of L ai p ie (n,T). 

There are in fact k + I free parameters in the definition of an element L(n, T) £ 



Tr a '2 e since the operator P(n, T) (see ( 3.30D ) only depends on the choice of the space 



Span{/W(n)}^j"o 1 (= KerP(n,T)), and not on the choice of the individual func- 
tions /W(n). Using again the linear independence of {(^(n)}^ 1 U {^±\n)}jZ} Q , 
the choice of span is equivalent to a choice of flags 

VqCViC . . . CVk-i and WqCWxC . . . C^h 

where V z = Span{/M(n)}; =0 and Wj = Span{/M(n)}^, cf. Q. 

The relations ( ]2.19| ) and ( |2.20 ) for L a jj :£ (n,T) imply similar relations for the 



sets 2? 

(3.36) 
(3.37) 



(fc,0 

a,/3;e 



Ad 



(-l) nI, /W 



= Ad 



(-i) 



= _p(*.0 



Here, in addition to (2.19) and (2.20), we use that the change of parameters a — > 
—a, /3 — » — /3, e — > e+a+ZJ exchanges tp^ (n) with tj)+\n) an d ^^(ti) with ip^}{ri). 
Analogously the change of parameters a —y (3, /3 — > a, e^e exchanges (n) with 
(-lyV^V) and ^(n) with (-l)"?/>-V)- 
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The Darboux maps between Jacobi functions (operators) represented by the first 
identities in (2.24), (2.25) and Remark 3.3 imply the following inclusion relations 



(3.38) 
(3.39) 



Ad: 



a,/3 — l;e a,/3;e" 

The function (n + e + a + /3)/(2n + 2e + a + (3) is the leading coefficient of the the 
operators D"(n,T) and D^_(n,T), see Sect. 2.3. Recall that the operator P(n,T) 
is normalized to have leading coefficient f . 



Remark 3.7. Note that fl3.3l| ), fl3 32| ) imply that for the operator P(n,T) ( |3.30[ ) 



defining an element L(n, T) in "D*£w. E the endomorphism L a ,p- £ (n, T) on KerP(n, T) 
has two Jordan blocks with eigenvalues f and — f and lengths k and Z, respectively. 
Insisting on multiple blocks with equal eigenvalues does not produce larger sets 
of transformations since the operator L a / j. £ (n,T) has a two dimensional kernel. 
Allowing k > \a\ or I > [3 in the cases a e Z or /? £ Z causes the operators P(n, T) 
and L(n,T) to have nonrational coefficients which does lead to bispectrality of 
L(n,T) as was noted in the introduction. 



At the end of this subsection we compute explicitly the coefficients of the oper- 
ators L[n, T) in V { ^ £ . Set 

(3.40) L(n,T) = a(n)T + b(n) + c{n)T- 1 

for some functions a(n), b(n), and c(n) (the dependence on A, B, C, D will not be 
shown). For convenience we denote the coefficients of the operator L a ^ e (n,T) by 
a (n), b (n), and c (n) : 

(3.41) L a ^ e {n, T) = a (n)T + b {n) + cqHT" 1 



(cf. eq. (2.14) for their values). Set also 

(3.42) det_ r (n):=det(/«(n+j))i=o,...,fc+l-i for r = 0, . . . ,k + l. 

j=—k—l,... ,0 

Note that 

(3.43) deto(n) = det(n) and det/ £+ /(n) = det (n + 1) = det(n + 1). 
Expanding the determinant ( 3.30| ) defining P(n,T) along the last column gives 



k+l 



(3.44) 



. det_ r (n) 
det(n) 



Proposition 3.8. The coefficients a(n), b(n), and c(n) of an operator L(n,T) £ 
^a'p-e are expressed in terms of the coefficients ao(n), bo(n), andco(n) of L a ^.^ e (n,T) 
and the functions functions f"> (n) (see eq. (3.42)) by the following formulas 

(3.45) a(n) — ao(n), 

ti ar\ u \ k i \a i N det-i(n + l) , ,det_i(ri) 

(3.46) b(n) — b (n) + a (n)— - — — a (n — 1)- 



(3.47) 



c(n) = co(n — k — I) 



det(n+l) ~' det(n) 

det(n - l)det(n + 1) 
(det(n)) 2 ' 
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Proof. We compare the coefficients of T and f in (3.33) and use eq. (3.44) for the 
operator P(n,T). This gives the formulas 

a(n) = ao(ra), 

det_i(n + l) , , det_i(n) 

o(n) — a{n)— — — — = b (n) — a (n — 1)- 



det(n + 1) 



det(n) 



which are equivalent to eqs. ( 3.45 ), ( 3.46 ). 

Similarly comparing the coefficients of j , - fe -'- 1 j n ( 3.33 ) gives 

c(n) dct (fc+l) (»-l) = Co{n k _ dBt_ (W) (») 



det(n - 1) 



det(n) 



which implies (3.47), taking into account (3.43). 



□ 



4. BISPECTRAL DARBOUX TRANSFORMATION AND AN INVOLUTION 

This section is a preparation for the next one where we show that the difference 
operators from are bispectral under some natural conditions on a and (3. 

Our proof is based on a result of || on Darboux transformations that preserve the 
bispectral property. Its application to the situation under consideration is nontrivial 
and requires an intrinsic characterization of a certain space of difference operators. 
This is done in terms of an involution of the algebra of difference operators with 
rational coefficients. 

4.1. A theorem on bispectral Darboux transformations. For a fixed choice 
of the parameters a, (3, e we define B a ^- e as the algebra of difference operators 
S(n, T) with rational coefficients for which there exists a differential operator G(z, d z ) 
(also with rational coefficients) such that 

(4.1) S(n, T)p a £ >P(n, z) = G(z, d z )p^(n, z). 

The set of all such operators G(z, d z ) is an algebra which will be denoted by B' a g. e . 
It is clear that 

(4.2) b(S(n,T)):=G(z,d z ) 
correctly defines a map 

(4-3) b : B.,, , - b',',., 

which is an antiisomorphism of algebras. In this setting eqs. ( [2.15] ), ( [2.16 ) mean 



that X e (n),L a ,ff ;E (n,T) g B a ,p;e, z,B a ^(z,d z ) g B' aJ3 . e , and 

(4.4) b(X e (n))=B ail3 (z,d z ), 

(4.5) b(L atf} . e (n,T)) = z. 

The triple (B a ,i3;e, B' a ^. £ , 6) is an example of a bispectral triple in the sense of 
Denote 

(4.6) K , fe = B a , fe nC(n), 

(4.7) ^ = ^ncw, 
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where C(n) and C(z) stand for the algebras of rational functions in the variables n 
and z, respectively. Let 

-4a : /3;e = 6 _1 {K,' a ^. e ) , 



(4.8) 
(4.9) 
It is obvious that 
(4.10) 
(4.11) 
and 
(4.12) 
(4.13) 



X 



,,. i;e - (k<>- 

= C[z], 



a,f3;e 



-4a,/3;e^C[L Qj/3;e (n, T)}. 



Later in Remark 4.3 we will show that the inclusions in ( 4.12 ) and ( 4.13 ) can be 
strengthen to give two equalities. 

As was noted in Sect. 3.1, if a difference operator q(L a ^^ e {n, T)) G A a ,f3; e S 
C[x]) is factorized as a product of two operators Q(n,T) and P(n,T) 

q(L a ^. e (n,T)) = Q(n,T)P(n,T), 

then the function 

V(n,z) = P(n, T)p^(n,z) 

is an eigenfunction of the difference operator P(n,T)Q(n,T) : 

P(n, T)Q(n, T)V(n, z) = q{z)^{n, z). 

We will give a version of Theorem 1.2 from || which provides general sufficient 
conditions on the operators P(n,T) and Q(n,T) under which ^(n, z) is also an 
eigenfunction of a differential operator in the variable z. (The original result of |j| 
deals with "bispectral" Darboux transformations in an arbitrary associative algebra 
but in the form to be used, needs an additional refinement.) 

Theorem 4.1. Assume that the operator q(L a ./3 ;£ (n, T)) G A a ^- £ is factorized as 

(4.14) q(L a ^ £ (n, T)) = (Q (n, T^n)- 1 )^)- 1 P (n, T)) 

for some difference operators Po(n,T),Qo(n,T) G B a ^- e and rational functions 
/j,(n),v(n) G /C Qj/ 3 ;e , such that the coefficients of the operators fi(n)~ 1 Po(n,T), 
Qo(n, T)i/(n) -1 are correctly defined for n G Z. Then the function 

(4.15) Hf(n,z) = (/i- 1 (n)P («,T))p«^(n,2) 
satisfies the relations 

(4.16) (^(n)- 1 P (n, T))(Q (n, rMn)- 1 )*^, z) = g(z)*(n, z), 

(4.17) 6(Po)(«, S*)6(Qo)(js, d t )q{z)-H(n, z) = v(n)v(n)V{n, z), 
i.e. it is bispectral. 

we do not assume that the rational functions p{n) 



4.1 



Note that in Theorem 
and j/(n) _1 are well defined for n G Z, but only that the "ratios" p(n)~ 1 Po(n,T) 
and Qo(n, T)^(n) -1 are. Because of this a small modification of the original proof 
from H is necessary. 
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First of all since the algebra B' a p. E has no zero divisors, eq. ( 4.14 ) implies (see 
EJ) 

(4.18) M(z,d z ) (bfi)(z,d z ) = (bP ){z,d,)q(*)~ 1 (PQo){.z,d z ). 
For all values of n for which fi(n) does not vanish we have 

*(n, z) = n{n)-\bPfi{z, d z )p^(n, z) 

and (4.17) holds, as a consequence of ( p~18| ). The validity of ( p~17| ) for all n G Z 
follows from the definition (4.15) of ^(n, z) and the fact that p" ,/3 (n, z) has an 
expansion in z around z = 1 with coefficients that are rational functions in n 
(recall (gig)). 

Returning to the sets TJ^'p-e °f Darboux transformations from the operators 
L a ,p;e{n, T), we need to find which of the operators P(n, T) from eq. ( 3.3C ) can be 
expressed in the form fj,(n)~ 1 Po(n, T) with fj,(n) and Po(n, T) as above. According 

to Theorem 4.1 the corresponding operators L(n,T) G 2^fj^ s will be bispectral 
with bispectral eigenfunction (3.34) (see also ( 4.1 5[ ) ) . For this we need an invariant 
description of the linear space of difference operators 

(4.19) TZa : /3- s — Span{/i(n)~ 1 S'(n, T) \S(n,T) G B Q>l a ;e , fi(n) G IC a ^- e , such that 

^((n) _1 S(n,T) is well defined forn G Z}. 

This will be obtained in the next subsection. Here we would like to note that the 
dual object - the linear space of differential operators 

(4.20) K,p-,e = Span^z)-^, d z ) \ G(z, d z ) G B' a ^ £ , g(z) G K.' a>p . E } 

is much easier to describe. It is just the space of differential operators with rational 
coefficients. This is a consequence of the fact that the commutator 

[B a ,p{z, d z ), z] = 2(z 2 - l)d z + ((a-l3) + {a + /3 + 2)z) 
is a first order differential operator that belongs to B' a g. e and z G K! a p. e (see eq. 
( 4.10| )). Unfortunately for our proof of the fact that the operators from D^'p. e are 
bispectral we need the space IZa.p-ei an d n °t the space lZ' a p. e . 

4.2. Description of TZ at p- e . Denote by A the abstract algebra of difference oper- 
ators M(n,T) with rational coefficients; that is the algebra over C, generated by 
rational functions in n, the shift operator T, and its inverse T , subject to the 
relation 

Th(n) = h(n + 1)T, for all rational functions h(n). 

Here we do not require that the coefficients of an operator M(n,T) in A be well 
defined for n G Z. More explicitly these coefficients could have poles at some neZ. 
The subspace of A consisting of operators having this extra regularity property will be 
denoted by A reg . We will identify the space of difference operators with rational coef- 
ficients acting on functions / : Z — > C with A reg . In particular, S Qi /3 ;£ c7?.a,/3; e cA res . 
Define an involution / in the algebra A acting on rational functions h(n) by 

(4.21) (Ih)(n) = h{-{n + 1s + a + [3 + l)) 
and on the shift operator T by 

I(T) = T -1 . 
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The involution / is correctly denned since 

I(T) (Ih)(n) = (Ih)(n + l)I(T). 
Denote the fixed points of I in A by A 1 : 



{M(n, T) G A | I(M(n,T)) = M(n, T)}. 

x (e + a + l) n 
<j>(n) - 



(4.22) £ 
Let 

(4.23) 

(e + l) n 

(cf. the definition (|2.13|) oip^{n,z) for a £ Z <0 ). 



Theorem 4.2. TTie space of difference operators TZ a ,/3;E defined in (4.19) is char- 
acterized by 

(4.24) TZ aJj , e = M m (A 1 n A rc «) , 

i.e. a/ter conjugation by 4>{n)~ l all operators from r R- a ,j3;e are I -invariant. 
Proof. Consider first the case a £ Z<o- Let 

(4.25) p^(n,z) = <f>(n)- 1 p^(n,z). 
The expression (2.13) implies 

(4.26) Ps ,f} {n, z) = F(-(-n + e), n + e + a + [3 + 1; a + 1; (1 - z)/2). 

Let B a ,/3;e, B' a ^. e , K. a ,P;e, ^' a .fi-^ Ra&e, an d K' a p. s , denote the B, JC and 72. objects 
associ ated w ith t he func tions p"'^(n,z) (see the beginning of Sect. 4.1 and eqs. 
([0]), (|47|) , ( jrj ), (|4.2d|) for the appropriate definitions). Obviously 

and A: Q:/3;e = JCa,^, fi£ j/3;E = ^ !/3 . £ , K! a<p . e = K,p-,e = B Lp-e- In this 

notation, the statement of the theorem is equivalent to 



(4.27) 



To prove that the l.h.s. of (4.27) is contained in the r.h.s., let us fix an oper- 
ator R(n,T) G ~R- a ,f3-£- There exists a difference operator S(n,T) G B a> p- E and a 
function j2(ri) G IC a> p- E such that R(n,T) = /i(n) _1 S(n,T). We will prove that all 
operators from B a ,p; S are /-invariant. This in particular shows that all functions 
from /C aj( a ;e CiB aj( 3 ;e are /-invariant and so are all operators from TZ a ,/3-e- 

If S(n,T) G B a ^- e , then there exists a differential operator G(z,d z ) for which 

(4.28) S(n, T)p^(n, z) = G(z, d z )p^(n, z). 

The fact that the hypergeometric function F(a,b; c; x) is symmetric with respect 
to a and b, and formula ( [4.26D for p™'P(n, z) imply 

(4.29) I (S(n, T)) p^(n, z) = G(z, d z )p^(n, z). 
Combining ( fl2S| ) and (fD|), we conclude that 

(S(n, T)-I (S(n, T)) )p^(n, z) = 0. 
This is only possible if 

l(s(n,T))=S(n,T). 
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The harder part of the proof of (4.27) is to show that any /-invariant difference 
operator from A reg belongs to TZ a ^- £ . It is sufficient to prove that for any R(n, T) £ 
A 7 there exists S(n,T) s Ba,p-,e and A( n ) G fca,0;e such that 

R(n,T) = fi(n)- 1 S(n,T). 

First let us write formulas ( 2.15|) and ( 2.1 6| ) in terms of p" ,/3 (n, z). Eq. ( 2.16 ) 
remains unchanged: 

(4.30) X e (n)p^(n, z) = B aJj {z, d z )p^(n, z), 
while eq. ( 2.15 ) becomes 

(4.31) L a , 0;e (n, T)p^(n, z) = zp^(n, z) 
with 

(4.32) L a , 0;e (n, T) = ^{n)- 1 L a ^, e {n, T)</>(n) 

2(n + e + a + l)(n + e + a + (3 + 1) 



+ 



(2n + 2e + a + (3 + l)(2n + 2e + a + /? + 2) 

« 2 -/3 2 

(2n + 2e + a + (3)(2n + 2e + a + f3 + 2) 

2{n + e){n + e + P) , 



T 



{2n + 2s + a + f3){2n + 2e + a + f3+l) 

The algebra A has a natural Z>o hltration where A^ consists of all operators from 
A with support [-d, d]. Denote Aj=A'nA d . 

We will prove that any difference operator R^(n,T) £ A^, can be decomposed 
as a sum 

(4.33) R^n, T) = ftn)- 1 ^, T) + R^n, T) 
where 

(4.34) S(n, T) £ B a , 0;s , /i(n) £ K a , ft£ , 

(4.35) tiKTleAi, 

Since Aq = C(A e (n)) (any /-invariant rational function in n is a rational function 
in A e (n)), by induction on d eq. Q4.33D implies that 

R^n,T)eK at/3 . e . 

A straightforward computation yields 

(4.36) ad Ae( „)T d = (A e (n) - A £ (n + d)) T d 

= -d(2n + 2e + a + /? + d + l)T d , 

and thus 

ad Ae(ri) (ad Ae( „ ) + l)L a ,p. s = 2(Id + I)((n + e + a + l)(n + e + a + /? + 1)T). 

So 

(4-37) (ad AE(n) (ad AE( „) + l)Z ai ^ E ) 

= 2 d (Id + /) f JJ(n + e + a + i)(n + £ + a + (3 + i)T d ^j + U d -i. 
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for some U d -i 6 A^-i- (Here we use the 7-invariance of L Qj( g ;e (n, T).) Denote for 
simplicity 

d 

c d (n) = 2 d Y[((n + e + a + i)(n + e + a + l3 + i)) 



i=l 



and let 



R d (n,T) = fl ?rt r ' < n )Mn) G C[n]. 



Using (pL36b and (p~37j) we obtain 
(4.38) 



1 a + + d + l \ / 1 
a d (-^ad Ae(n) £ U d [-ad Ae(n) 



a + (3-d+l 



-ex 



, 1 , a + P-d+1 
x c d I ^ad Ae („) e 



)( 



2cT 2 
(Id + /) f & d ( n ) c rf (n) (Jc d )(n) 



adA e (n)(ad Ae („) + l)L a ,p. t 



a d{n) Td 



bd(n) 

for some other Ud-i G ^ I d _ 1 . There exists a polynomial q d (n) for which 

b d (n) {Ib d )(n) c d (n) {Ic d ){n) = q d (X £ (n)) 
because the polynomial in the l.h.s. is clearly /-invariant. Denote by S(n,T) the 



difference operator in ( 4.38) ). The l.h.s. of (1.3S) implies that S(n,T) belongs to 
Ba,0;s and the r.h.s. implies 

^(n,T)~(< Z (A £ (n))- 1 §(n,T)eA^ 1 



which completes the proof of Theorem 4.2 



□ 



Remark 4.3. Any fuction Jl(n) £ IC a ^- e is /-invariant and therefore is a rational 
function in A s (n). In fact Jl(n) should be a polynomial in X e (n). Indeed if = 
p(X e (n))/q(X e (n)) for two polynomials p (x), q (x) G C[x] such that q(x) J(p(x), then 
there exists a differential operator G(z, d z ) with rational coefficients such that 

G(z > d z )p^(n,z) = ^^p^(n,z) 

which implies 

(4.39) p{B a ,p{z, d z )) = G(z, d z )q(B a A z > 9,)). 

This is impossible; if zq is a root of q(x) and p(x) of multiplicities d\ > c?2, then 
there exist a holomorphic function g(z) in a domain f2cC such that 

(B a>/3 (z, d z ) - z Q ) dl g(z) = and q(B at p(z, d z ) - z Q )g(z) 



which contradicts with (4.39). Since /C aj( g ;e = IC a ^- e we finally obtain 

= C[A e (n)], 

^ i/3;£ = c[s Q ^(z,a z )], 



as promised following (4.12) and (4.13) 
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Remark 4.4. The second order bispectral differential operators of the even case 
of Duistermaat-Grunbaum's classification |6| are obtained as Darboux transforma- 
tions from the Bessel operators 



, 2 Kk-l) 1 



L k (x,d x ) = d x — , k e Z- 



in the sense of (3.1). More precisely for each operator L(x,d x ) of this family there 



exists a differential operator with rational coefficients P(x, d x ) such that 

L(x, d x )P(x, d x ) = P(x, d x )L k (x, d x ). 

In addition, the operator P(x,d x ) satisfies 

(4.40) P(x,d x )=P(-x,-d x ). 

Let X denote the involution of the algebra of differential operators with ratio- 
nal coefficients induced by the diffcomorphism x i— > — x of C (i.e. (XS)(x,d x ) = 



S(—x,—d x )). Then (4.4C) means that P(x,d x ) is invariant under X. This gives the 
relation of the approach of this paper via the involution / and the space TZ a ,p-e to 
the construction of Q . 

5. BlSPECTRALITY OF Z^ffl-e 

In this section we prove our main result: when the parameters a and [3 are subject 
to certain natural integrality conditions, the difference operators from T>^'g. e are 

bispectral. As an example, for each L(n, T) S £>2 o°e we ^ nc ^ a dual differential 
operator of order 10. 



The conditions (2.18) on a, /3, e are assumed throughout this section. 



5.1. Proof of the main result. The conjugation by the function 4>{n) (see (4.23)), 
used in Theorem [h^, leads us to consider the functions := <jJ± (n) / 4>(n), := 
^± (n) / 4>{n) . Because of eqs. ( 3.18 )— ( 3.21 ) they are explicitly given by the formulas 



(5.1) ^(n)^ ( ' (n + £)Mn + £ + ". + /j + 1) \ 

(5 2) * (4) (n)- (e + + l)n(e + l)n (-(n + e + a + /9))i(n + e + l) t 



(5.3) $ W (n) 



(e + a + l)„(e + a + /3 + l)„ (-a + l) 2 (-2) 1 

(g + /?+!)„ (-(n + £)). t (n + e + a + /j + l) j 
(-!)"(£ + a+l)„ (/3 + l)i2' 



rt - 4 n ^(i), s = (£ + !)» (-(n + £ + a + /3)) t (n + e+ 1), 

1 j (-l)«(£ + « + /3 + l)„ (-/3+l) t 2< 

Lemma 5.1. If a & Z, then for i < |a| — 1, $^(n) and (n) are I-invariant 
rational functions of n. 

If a 6 Z and/3 e Z then fori < \a\-l, j < <f>+(n), *+ (n), (-1)"$ 0) (n), 

and (—l) n ^!_(n) are rational functions ofn, $^(n), 'J/j! (n) are I-invariant, and 

(5.5) / ((-l)"$ 0) (n)) = (-l) a+/3 ((-l)"$ 0) (n)) , 

(5.6) / ((-l) n *^(n)) = (-1) Q+/3 ((-l)™* '^) . 



24 



F. A. GRUNBAUM AND M. YAKIMOV 



Proof. First note that 

(-(n + e))i(n + e + a + f3 + 1); 



i-l 



]f[(-(n + e) + r)(n + £ + a + /3+l + r) 

r=0 

i-l 

(-l) fc JJ(A(n) - r(a + /3 + 1 + r)) 



i-l 



and similarly 

(-(n + e + a + /3)),(n + e + 1) 4 = JJ(A(n) - (a + /? - r)(r + 1)) 

r=0 

are /-invariant polynomials in n. 

To prove the first statement of the lemma we use a similar computation. Re- 
stricting to the case a € Z>o : 

(£ + /?+ l)n(£+ l) n 

(e + a + l)„(e + a + /3+l)„ 
(e+l) a (e + /3 + l)a 



(n + e + l)a(n + e + /3 + l)a 

= (g+ l)q(g + /?+ l)q 

+ £ + r )( n + £ + /? + a + l-T) 

= (£+l) a (£ + /3+l) a 

The proof of the second statement is analogous. Assuming a, [3 e Z >0 and /3 > a 
we obtain 

(g + /3+ l) n = (g + a- + n + l)ff_ a 
(g + a + l)„ (e + a + l^-a 

_ Q(n) ]lr=i (« + e + a + r)(n + e + /5 + l-r) 



(e + a + 1)/3-q. 
<?(«) Jiff 1 (A(n) + (a + r)(/3 + 1 - r)) 



with 



(e + a+l)^ 

f 1 

g(n) 



1 , if (3 + a is even 

n + e + (a + /3 + l)/2, if /? + a is odd 

(Since a e Z, the first condition is equ ivalent to 2|(/3 — a) and the second one to 
2 /(/? — a).) To finish the proof of (5.5) we just observe that 

(5.7) I(n + e + (a + (3 + l)/2) - -(n + g + (a + + l)/2). 

The remaining cases for a, /3 £ Z are treated analogously. 
The identity (5.6) follows from the analogous formula 

(g+ l)n _ (g + 

(e + a + /3 + 1)„ ?(n) (A(n) + r(a + ^ + ! _ r) ) 
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and eq. (| 

Throughout this proof, for a real number x by [x] we denote its integer part. □ 

Theorem 5.2. Assuming ( |2.18 ), the following sets consist of bispectral difference 
operators: 

1.) D^afe if a EZ and k< \a\, 

2 -) v aie ifPeZandl<\P\, 

3.) X>£f if a, (3 £ Z and k < \a\, I < \[3\. 

When the conditions (2.18) are not met but the operator L a ^ ;e (n, T) is still well 
defined the arguments below can be adapted properly. We do not pursue that here. 



Proof. Because of the relation ( 3.37 ) the second case follows from the first one. 

Let us restrict to instances 1 and 3 of the theorem above. In each of them we can 
assume that k + l is even using ( |3.38 ). Fix an operator L(n, z) e T^'b-ei determined 
by a choice of the functions {/^ (ft) i.e. a choice of admissible values of the 

complex parameters A,B,C,D (see Sect. 3.3). It has the eigenfunction ^(n,z) 
defined in ( |3~3^ ) 

L(n, T)^(n, z) = z^(n, z), 

cf. ( 3.35 ). We need to show that there exists a differential operator B(z, d z ) having 
^(n, z) as an eigenfunction, that is 

B(z, d z )V(n, z) = 6(n)V(n, z) 

for some function 9{n). 
Define the functions 

F^(n) = f (l) (n)/(j>(n), i = 0, . . . , k + I - 1. 

Let us put s := (k + l)/2 and consider the operator 

i^°)(n-s) ... F( k+l -V(n- 



(5.8) 



P(n,T) = (-l) 



■nl 



T~ 



F(°)(n + s) 



F( k+l -V(n + s) T s 



It is a regular difference operator with kernel given by Span{F«(n)}*lt"o _1 . Hence 
it is related to the operator P(n,T) (recall eq. ( 3.30| )) by 



(5.9) 
where 



P{n, T) = d{n)- 1 (f){n)- 1 T- s P{n,T)(j){n) 



d(n) = (-l) ni det(F«(n + j))^" 1 '" 1 



)i,j=0,-k-l 

is the leading coefficient of T~ s P(n,T). Lemma [5.l| implies that F^(n), ... , 
ir(*-i)( n ) 5 and (-l) n F ( fc )(n), ... , (— l) n F(* + * -1 )(ra) are rational functions in n. 
This implies that for i = k, . . . , k + / — 1 and for all j G Z, (— l) n F^(n + j) are 
also rational functions in n and thus P{n, T) has rational coefficients. In addition 
Lemma 5.1 gives 



/(>«(n + j)) =F^(n-j),i = 0,. 



,fc-l,jeZ 



and 



((-1)"F« (n + j)) = (-I)(^)/ 2 ((-1)"F« (n - j)) , i = k, . . . , k + I - 1, j e Z. 
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Taking into account that I(T) = T 1 we obtain 



(5.10) 



I[P(n,T) 



where the factor (— l) s comes from exchanging the pairs of rows (l,k + l+ 1), 
(s,s + 2). Set 



q(ri) 



1, if s + (a + /3)Z is even 

(n + £ + (a + /? + l)/2), if s + (a + /3)Z is odd 



and consider the operator 
(5.11) 



P(n,T) = q(n)P(n,T). 



Because of the conditions ( [2.18 ), q(n) does not vanish for n G Z. Taking into account 
(^9|) one sees that P(n, T) is related to P(n, T) by 



(5.12) 



P(n,T) = y^-fi r-^w-^r)^). 



g(n - s)0(n) ' 

Since T) is a regular difference operator and 4>(n) does not vanish for n G Z 
(recall ( 2.1 8| ) ) , there exists a difference operator with rational coefficients Q(n,T) 
such that 

(L a ^ £ (n,T) - l) k (L a ^, £ (n,T) + l) 1 = (0(n)- J Q(n, T)0(n)) (^n)" 1 ^, T)0(n)) 



From eqs. (5.1C) and (5.11) it follows that P(n, T) is /-invariant. Finally combin- 
ing this with the J-invariance of <f>(n)~ x L a> p. e (n, T)<p(n) = L a> p- e (n,T) (see (4.32)) 
implies the i-invariance of the operator Q(n,T). Theorem 4.2 now gives 

(?!>(n) _1 P(rz.,T)0(7i), ^(n)" 1 ^, T)^{n) G TZ a ^, £ . 

Applying Theorem [4.l| we obtain that the function 

(5.13) tf (n, z) = <f>(n)P{n, T)<t>{ri)~ r p%' p (n, z) 
is an eigenfunction of a differential operator S(z, d z ) 

(5.14) 5(z, 5 a )*(n, z) = h(\ E (n)y$(n, z), 



for some polynomial h(x). Because of (5.12) our original function ^(n, z) G "D^ i} £ 
is related to ^(n,z) by 



(5.15) 



*(n,z) = P(n,T)p^(n,z) = d [ n ^ n Jl j-^/^fn, z). 

q(n - 8)<p(n) 



Eq. (5.14) implies that ^(n, z) is an eigenfunction of the same operator B(z,d z ) 
with eigenvalue T- < - k+l ^ 2 h(X £ (n)) : 



(5.16) 



B(z, <9 z )*(n, z) = h(X £ (n -(k + J)/2))*(n, z) 



□ 
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5.2. An example: the set l^oe- -"- n this final subsection we consider in detail 
the case a = 2,/3 = 0,fc = 2,Z = and use this example for two different 
purposes. First we give the reader a guided tour through the results in this paper: 
we start with the function p^'°(n, z) f rom J2~lg| ), give the ingredients needed to 
build th e diff erence operator P(n, T) ( |3.30| ) and the corresponding eigenfunction 
^{n, z) ( |3.34 ), and end with a description of the strategy used in the construction 
of a differential operator in the variable z giving a bispectral situation. The algebra 
of possible differential operators in z contains some whose order is lower than the 
one resulting from this construction. We close this subsection with an explicit 
expression for the (essentially unique) bispectral operator of minimal order and 
material related to this operator. 



The functions ipl (n) and ip+\n) (i = 0, 1) from ( |3.1<j )-( 3.19 ) are given by 



(5.17) 
(5.18) 
where 



^(n)= (n + g K +1)a , 
<>(n) 



{n + e + 1) 2 ' 



2" 



(5.19) «=(e+l)(e + 2). 

The conditions ( [2.18; ) reduce to e £ Z. 



An element L{n,T) € T>2q\ is determined by a choice of the functions 



f^(n)=A ^>(n)+B ^>(n), 



,(0) 



(0), 



/W(n) 



Brff (n) + A ^(n) + B ^{n), 

cf. Sect. 3.3. We will restrict to the generic case when Ao 7^ 0. In this case 
we can assume that Aq = 1 and Ai = by dividing f^°'(n) by and then 
subtracting from f^ 1 \n) the term Ai/W(n). Recall that L{n,T) depends only on 
Span{/(°) (n), (n)}. Once this space has been spec ified by the choice of Bo, B\ we 
can build the difference opera tor P(n,T) as in ( 3.3C ) and we get the eigenfunction 
*(»,z) of L(n,T) from (fljj). 

The theory developed in Sections 4 and 5 makes it convenient to introduce the 
difference ope rator s P(n,T), see (E]l), and P{n,T), see ( 5.11 ), related to P(n,T) 
by (HJ) and (EH)- 

The m ain point i n th e proof of Theorem |5.2| is that the operator P(n, T) defined 
in ( 5.11 ) (see also ( |5.8| )) is /-invariant and thus 4>(n)P(n,T)(f>(n)^ 1 e 72.2,0;e- This 
implies that the function 

*(n,«) = <?!)(n)P(n,r)0(n)" 1 p"' /3 (w,2) 
(see J5T^) ) can be expressed as 

(5.20) ¥(n, 2) = ^(n) _1 G(js, d,)p? fi (n, z) 

for some differential operator with r ational coefficients G(z,d z ) and some polyno- 
mial /i(n) (recall the definition (4.19) of 72-2,0;e)- Now any operator B(z, d z ) that is a 
Darboux transformation from h(B2fi{z, d z )) for some h(x) £ C[x] via the operator 
G(z,d z ), i.e. 



-,(!)/ 



(5.21) 



B(z, d z )G(z, d z ) = G(z, d z )h(B 2 , (z, d z )) 
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will satisfy 



B(z,d z )V(n,z) = h(\ e (n))V(n,z) 



(a differential analog of ( |3.3[ )). The function ^(n, z) is related to ^{n, z) by ( |5.15 ) 
and is also an eigenfunction of B(z,d z ) but with eigenvalue h(\(n — 1)) 

B(z, d z )^(n, z) = h(X e (n - l))*(n, z), 

see ( |5.16| ). Combined with ( |3.35| ) 

L{n,T)^(n,z) = z^(n,z) 



this gives the desired bispectral pair (L(n, T), B(z, d z )). 

The J-invariance of the operator P(n, T) in this special case can be observed 
directly. Because of (|5.17|) , fl5.1S|) the functions F^(n) = f^(n)/<p(n), i = 0,1, 
see ( |4.23 ), are given in terms of 



by 



F«(n) 



= 1 



A E (n) = (n + e){n + e + 3) 



BoA £ (n) 



Gk 

BiAe(n) , 



6k (^e( n ) 

The operator P[n, T) is given by 



(A £ (n) 



So 

2 



P(n,T) - (n + e + 3/2) 



F<-°\n-l) F«(n-1) T" 1 
F(°)(n) 1 



and it is /-invariant because of the /-invariance of A e (n) and the skew invariance of 
the factor in front compensating the effect of the exchange of first and third row. 



An operator G(z, d z ) satisfying (5.20) is generated from the proof of Theorem 4.2 
It is of high order and the one of minimal order 10 has the following form 

G{z, 8 Z ) = {z- lf{z + Ifd™ + (z- l) 5 (z + l) 4 (57z + 7)d 9 z + 

+ A{z - l) 4 (z + l) 3 (311z 2 + 68z - 43)5f + {W q k 2 (z - l) 2 (z + l) 2 + 
+ 2(18793z 4 + 5796z 3 - 15734z 2 - 3636z + 1501))dJ + . . . 



Theorem 4T guarantees that ( 5.21 ) is satisfied for some polynomial h(x). It also 
generates such a polynomial but it is again of high order. The one of minimal order 
5 is given by 

h(x - 2) = x 5 - 5x 4 + (10B Q K 2 + 8)x 3 - 

- (30Bik 2 + 20B o k 2 + 4)a; 2 - IhBl^x. 
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Given G(z,d z ) eq. ( 5.21 ) determines the dual bispectral operator B(z,d z ) of 
L(n, T) of minimal order uniquely. It is given by 

B(z, 8 Z ) = {z- lf{z + lfdf + 50(z - l) 4 z(z + lfdl + 5{z ~ l) 3 (z + l) 3 x 
x (llz - 5)(17z + 7)<9f + 160(z - l) 2 (z + l) 2 (52z 3 - 7z 2 - 28z + l)d 7 z + 
+ (30B 2 k 4 z + 120B 1 k 2 (z - 1) + 120B k 2 )^ + (180B k 2 (z - l) 2 z(z + l) 2 + 
+ 240(z - l) 2 (337z 3 + 504z 2 + 141z - 30))df + (-30B 1 k 2 (z - l) 2 (z + l) 2 + 
+ 120S o k 2 (z - l)(z + l)(8z 2 - z - 3) + 120(z - l) 2 (641z 2 + 758z + 161))<9 4 + 
+ (-240 j B 1 k 2 (z - l)z{z + 1) + 240S k 2 (7z 3 - 3z 2 -7z + l) + 960(z - l) 2 x 
x (26z + 19))9f + (-60B 1 k 2 (z - l)(7z + 5) + 120B k 2 (2z + l)(3z - 5) + 
+ 1440(z - lf )d 2 z - (30B^k 4 z + 120Bik 2 (z - 1) + l2QB n 2 )d z . 

In the cases k = 1, I = 0, 1 and £ = the dual bispectral operator of minimal 
order was determined in [|(| . 



[l 

[2. 

[3: 
p: 

p: 

[9 
[10 

[11 

[12; 

[13; 

[14: 
[15; 

[is: 
[17: 
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